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Abstract. I discuss recent progress in calculating quarkonia correlators and spectral 
functions on the lattice in relation with the problem of quarkonia dissolution at high 
temperatures and heavy quark transport in Quark Gluon Plasma. 
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1. Introduction 

Heavy quarkonia play an important role in studying hot and dense strongly interacting 
matter. Because of the heavy quark mass quarkonia binding can be understood in 
terms of the static potential. General considerations suggest that quarkonia could melt 
at temperatures above the deconfinement temperature as a result of modification of 
inter-quark forces (color screening). It has been conjectured by that melting of different 
quarkonia states due to color screening can signal Quark Gluon Plasma formation in 
heavy ion collisions pQ . Many studies of quarkonia dissolution rely heavily on potential 
models [3 El IH ISl El IZ| • However it is very unclear if such models are valid at finite 
temperature jH]. 

The problem of quarkonium dissolution can be studied more rigorously in terms 
of meson (quarkonium) spectral functions. Lattice calculation of charmonium spectral 
functions appeared recently and suggested, contrary to potential models, that J/ip and 
Tjc survive at temperatures as high as 1.6Tc jHl CHI HH 1121 • ^as been also found that Xc 
melts at temperature of about l.lTc fUElEl- There are also preliminary calculations 
of the bottomonium spectral functions [T^ I15j. 



2. Meson correlators and spectral functions 

In lattice QCD we calculate correlators of point meson operators of the form 

J//(t, x) = q(t, x)THq(t, x), (1) 

where Th = 1, 75, 7;^, 757^, 7/^7;/ and fixes the quantum number of the channel to scalar, 
pseudo-scalar, vector, axial-vector and tensor channels correspondingly. The relation of 
these quantum number channels to different meson states is given in Tab. ^ 
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Table 1. Meson states in different channels 



Most dynamic properties of a finite temperature system are incorporated in the 
spectral function. The spectral function (Jh{Po-,P) for a given mesonic channel if in a 
system at temperature T can be defined through the Fourier transform of the real time 
two-point functions and D'^ or, equivalently, as the imaginary part of the Fourier 
transformed retarded correlation function [16], 

^h{Po.P) = ^iD>ipo,p) - D<{po,if)) = (2) 
-ImD§{po,p) 

IT 

Dp^\po,P)= I ^,e^'-D>J<\xo,x) (3) 
D^{xo,x) = (Jh(xo,x) Jh(0,0)) 

D<{xo,x) = {Jh{0,6) Jh{xo,x)), Xo>0 (4) 
The Euclidean time correlator calculated on the lattice 

GHir,i>) = J £xe'P^-''{T,MT,x)MO,0)) (5) 

is an analytic continuation of the real time correlator Gh{t,p) = D^{—iT,p). 

Using this equation and the Kubo-Martin-Schwinger (KMS) condition for the 
correlators 

D^{xo,x) = D<{xo + i/T,x), (6) 

one can relate the Euclidean propagator Guij^p) to the spectral function in Eq. 0, 
through the integral representation 

/"OO 

G{t,p) = / duja{uj,p)K{ijj,T) 
Jo 

cosh(^(r-l/2T)) 
sinh(c./2r) • 

To reconstruct the spectral function from the lattice correlator G{t, T) this integral 
representation should be inverted. Since the number of data points is less than the 
number of degrees of freedom (which is (9(100) for reasonable discretization of the 
integral ) spectral functions can be reconstructed only using the Maximum Entropy 
Method (MEM) J7j. In this method one looks for a spectral function which maximizes 
the conditional probability P[(j\DH] of having the spectral function a given the data D 
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and some prior knowledge H which for positive definite spectral function can be written 

as 



P[a\DH] = exp( 



where 
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duj 



cr(co') — m{(jj) — cr(co') ln( 



m[u] 



(8) 



(9) 



is the Shannon - Janes entropy. The real function m{uj) is called the default model 
and parametrizes all additional prior knowledge about the spectral functions, such as 
the asymptotic behavior at high energy In order to have sufficient number of data 
points either very fine isotropic lattices [HI [T21 US] or anisotropic lattices [01 Ell [13 [H] 
have been used. 



3. Charmonia correlators and spectral functions 

The spectral function for pseudo-scalar charmonium spectral functions calculated on 
anisotropic lattice is shown in Fig. 1. The first peak in the spectral function 
corresponds to rj^ilS) state. The position of the peak and the corresponding amplitude 
(i.e. the area under the peak) are in good agreement with the results of simple 
exponential fit. The second peak in the spectral function is most likely the combination 
of several excited states as its position and amplitude is higher than what one would 
expect for pure 2S state. The spectral function becomes sensitive to the effects of the 
finite lattice spacings for uo > 5GeV. In this u region the spectral functions becomes also 
sensitive to the choice of the default model. This is because only a very few data points 
in the correlator carry information about the spectral function in the region u > 5GeV. 

Also shown in Fig.l is the spectral function in the scalar channel from Ref. ^31- The 
1st peak corresponds to Xco(l-P) state. The correlator is more noisy in the scalar channel 
than in the pseudo-scalar one. As the results the Xco(l-P) peak is less pronounced and 
has larger statistical errors. The peak position and the area under the peak is consistent 
with the simple exponential fit. As in the pseudo-scalar case individual excited states 
are not resolved and the spectral function depends on the lattice spacing and default 
model for uj > 5GeV. Similar results have been found for the vector and axial-vector 
channels which correspond to J/ip and Xa states respectively. 

We would like to know what happens to different charmonia states at temperatures 
above the deconfinement temperature Tc. With increasing temperature it becomes more 
and more difficult to reconstruct the spectral functions as both the number of available 
data points as well as the physical extent of the time direction (which is 1/T) decreases. 
Therefore it is useful to study the temperature dependence of charmonia correlators 
first. From Eq. ([7|) it is clear that the temperature dependence of charmonia correlators 
come from two sources: the temperature dependence of the spectral function and the 
temperature dependence of the integration kernel K{t, u, T). To separate out the trivial 
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Figure 1. Charmonium spectral function in the pseudo-scalar channel (left) and the 
scalar channel (right) at different lattice spacings and zero temperature from Ref. 
m{u!) = 1 was used as the default model. 
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Figure 2. The ratio G{T,T)/Grecon{T,T) of charmonium for pseudoscalar channel at 
at a^'^ = 14.11GeV (left) and scalar channel at at a^^ = 8.18GeV (right) at different 
temperatures 



temperature dependence due to the integration kernel, following Ref. [12] for each 
temperature we calculate the so-called reconstructed correlator 

POD 

Greconir,T)= duja{uj , T = 0)K {t, UJ , T) . (10) 
Jo 

Now if we assume that there is no temperature dependence in the spectral function, 
then the ratio of the original and the reconstructed correlator should be close to one, 
G{r,T)/Grecon{'T,T) ~ 1. This way we can identify the cases when spectral function 
itself changes dramatically with temperature. This gives reliable information about 
the fate of charmonia states above deconfinement. In Fig. |21 we show this ratio for 
pseudo-scalar and scalar channels correspondingly calculated on anisotropic lattice [13]. 
From the figures one can see that the pseudo-scalar correlators shows only very small 
changes till 1.5Tc indicating that the rjc state survives till this temperature with little 
modification of its properties. On the other hand the scalar correlator shows large 
changes already at lAGT^ suggesting strong modification or dissolution of the Xco state 
at this temperature. 
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Figure 3. Charmonia spectral function in tlie psuedoscalar channel at a^^ = 
14.11GeV (left) and the scalar channel (right) at a^^ = 8.18GeV for zero and finite 
temperature [rgj . In the pseudo-scalar channel m(w) = 0.1 has been used as the default 
model. In the scalar channel we have used m{uj) — 0.038^^ as the default model as 
well as m{uj) = 1 for T LlGTc- 



More detailed information on different charmonia states at finite temperature can be 
obtained by calculating spectral functions using MEM. The results of these calculations 
are show in Figs. El Because at high temperature the temporal extent and the number of 
data points where the correlators are calculated become smaller the spectral functions 
reconstructed using MEM are less reliable. To take into account possible systematic 
effects when studying the temperature modifications of the spectral functions we 
compare the finite temperature spectral functions against the zero temperature spectral 
functions obtained from the correlator using the same time interval and number of data 
points as available at finite temperature. We see that spectral function in the pseudo- 
scalar channel show no temperature dependence within the statistical errors. This is 
in accord with the analysis of the correlation functions. Also the spectral functions 
show very little dependence on the default model. Similar conclusion has been made in 
Ref. fn\ IT^ where correlators and spectral functions have been calculated on very fine 
isotropic lattices as well as in Ref. jTU] where anisotropic lattice have been used. The 
pseudo-scalar spectral function was found to be temperature independent also in Ref. 
|j9j where correlators of extended meson operators have been studied on anisotropic 
lattices. The study of the charmonium correlators with different spatial boundary 
conditions provides further evidence for survival of the IS charmonia states well above 
the deconfinement transition temperature [T8] . 

The scalar spectral function shows large changes at 1.1 6Tc which is consistent with 
correlator-based analysis. Also default model dependence of the scalar correlator is large 
above the deconfinement transition (c.f. Fig. 3, right). This means that the Xco i^Po) 
dissolves at this temperature. Similar results for the scalar spectral function have been 
reported in fn\ IT^ . The results for the axial correlators and spectral functions are 
similar to scalar ones [12] as expected. 

The vector correlator, however, has temperature dependence different from that 
of the pseudo-scalar channel [TH]. This is due to the fact that the vector current is 



Quarkonia correlators and spectral functions from lattice QCD. 



6 



conserved and there is a contribution to spectral functions at very small energy c<j ~ 
corresponding to heavy quark transport fl^ 1^ . The transport peak in the spectral 
functions can be written as pUj 



where r] = T/{MD) with D being the heavy quark diffusion constant. Furthermore 
Xs{T) is the charm or beauty susceptibility and M is the heavy quark mass. To the 
first approximation the transport contribution to the spectral function gives rise to a 
positive constant contribution to the correlator Giow{t) — XsiT)T'^ / M j^HI resulting in 
an enhancement of the finite temperature correlators relative to the zero temperature 
ones, in agreement with the lattice data presented in P^ l^. Finite value of the diffusion 
constant D will give rise to some curvature in Giow{'t)- The smaller the value of D is, 
the larger is the curvature in Giowij)- Thus extracting Giowij) from lattice data and 
estimating its curvature can give an estimate for D. This, however, requires very precise 
lattice data which are not yet available |20j|. 

4. Bottomonium correlators and spectral functions 

Bottomonium correlators and spectral functions have also been studied using anisotropic 
[Ti] as well as very fine isotropic lattice with lattice spacing = 9.72GeV 
These studies, however, are far less detailed than the charmonium studies presented 
above and are still preliminary. In Fig. H] the temperature dependence of the ratio 
G/Grecon is showu for pseudo-scalar channel at different temperatures calculated using 
anisotropic lattices. This ratio shows almost no temperature dependence till 2.3Tc. 
This is expected as the pseudo-scalar 15* bottomonium state, the rjb, is much smaller 
than the corresponding charmonium state, thus survives till much higher temperatures. 
Also shown in Fig. |3]is the pseudo-scalar function at different temperatures. The first 
peak corresponds to the r]h state and survives above the deconfinement temperature in 
agreement with the analysis of the correlation functions. Other peaks are artifacts of 
the finite lattice spacing and MEM analysis. 

An interesting question is what happens to IP bottomonia states. They have sizes 
similar to the IS charmonia states and thus are expected to survive in the deconfined 
phase till temperatures of about 1.5Tc. In Fig. Ell show the temperature dependence of 
G/Grecon in the scalar channel corresponding to Xbo- As one can see from the figure the 
scalar correlator shows dramatic change across the deconfinement temperature and its 
behavior is similar to the behavior of the scalar correlator in the charmonium case. Also 
shown in Fig. Elis the bottomonium spectral functions in the scalar channel. Contrary to 
the pseudo-scalar channel the scalar spectral function shows significant changes above 
the deconfinement temperature and therefore it seems that the Xb state is strongly 
modified or even dissolved at 1.15Tc < T < 1.5Tc. Lattice calculations presented in 
Ref. jEj show similar results. In particular, they also show large increase of the scalar 
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Figure 4. Bottomonia correlators (left) spectral functions (right) in pseudo-scalar 
channel for different temperatures from Ref. |14| . 




Figure 5. Bottomonia correlators (left) spectral functions (right) in scalar channel 
for different temperatures from Ref. |14j . 



bottomonium correlator and strong modification of tlie corresponding spectral function. 



5. Conclusions 



In this contribution it has been shown how lattice calculation can provide information on 
quarkonia properties at finite temperature. The temperature dependence of the pseudo- 
scalar correlators as well as the spectral functions extracted using MEM shows that the 
15* charmonia states exist as a resonance in the deconfined phase till temperatures as 
high as 1.5Tc. On the other hand lattice calculation show that the IP charmonia states 
dissolve at T ~ I-IT^. Study of the bottomonia at finite temperature is also in progress 
[TH IT5] . It should be stressed that all lattice calculation discussed so far have been done 
in quenched approximation, i.e. neglecting the effect of sea quarks. To make contact 
with heavy ion experiments certainly the effect of the sea quarks has to be included, but 
computationally this is very expensive. Recent attempts to study charmonia properties 
at finite temperature in full QCD (i.e. with sea quarks) are reported in Ref. [2Sj. The 
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findings of Ref. 23 are consistent with tlie quenclied results. 

Recent studies of quarkonium properties at finite temperature using potential model 
claim agreement with the lattice data [3 El IZ| • The dissociation temperature of different 
quarkonia states quoted in these studies and defined as temperature where binding 
energy becomes zero is indeed significantly higher than before. However, potential 
models with screened potential also predict modification of quarkonia properties, which 
in turn lead to changes in the spectral function and correlators. Such changes are not 
observed in the lattice correlator fH] El] . Thus it is not clear at the moment whether 
modification of quarkonia properties at finite temperature can be understood in terms 
of a screened temperature dependent heavy quark potential. 
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